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ABSTRACT 

Aims. This work aims to provide a theoretical formulation of Surface Brightness Fluctuations (SBF) in the framework of probabilistic 
synthesis models, in which there are no deterministic relations between the different stellar components of a population but only 
relations on average, and to distinguish between the different distributions involved in the SBF definition. 

Methods. By applying the probabilistic theory of stellar population synthesis models, we estimate the shape (mean, variance, skew- 
ness, and kurtosis) of the distribution of fluctuations across resolution elements, and examine the implications for SBF determination, 
definition and application. 

Results. We propose three definitions of SBF: (i) stellar population SBF, which can be computed from synthesis models and provide 
an intrinsic metric of fit for stellar population studies; (ii) theoretical SBF, which include the stellar population SBF plus an additional 
term that takes into account the distribution of the number of stars per resolution element tff(N); theoretical SBF coincide with Tonry 
& Schneider (1998) definition in the very particular case that tf/(N) is assumed to be a Poisson distribution. However, the Poisson 
contribution to theoretical SBF is around 0.1% of the contribution due to the stellar population SBF, so there is no justification to 
include any reference to Poisson statistics in the SBF definition; (iii) observational SBF, which are those obtained in observations that 
are distributed around the theoretical SBF. Finally, we show alternative ways to compute SBF and extend the application of stellar 
population SBF to defining a metric of fitting for standard stellar population studies. 

Conclusions. We demostrate that SBF are observational evidence of a probabilistic paradigm in population synthesis, where inte- 
grated luminosities have an intrinsic distributed nature, and they rule out the commonly assumed deterministic paradigm of stellar 
population modeling. 

Key words. Galaxies: star clusters - Galaxies: stellar content - Methods: data analysis 

1. Introduction By construction, EPS needs to assume the existence of a 

probability distribution function that describes the mean value 

. Stellar population synthesis studies aim to decompose the in- relative contribution of different populations to the inte- 

tegrated light, Ltot;.i(A^*), of a sample of N* stars into a com- grated light, fi[(ni)/N*. In addition, it assumes that each popu- 

H bination of particular stellar classes or populations, «,-. These j^tjon ^an be characterized by a mean luminosity, ju'. (/,.,,)• EPS 

. .(^ studies can be developed by following two different methods, studies therefore refer to the generic emission of the integrated 

The first method ( jSpinrad & Tayloij | 1971|; |Fabeij | 1972t as ex- jjght of a set of systems rather than the partic ular emi ssion of 

amples) completes a direct decomposition of the integrated light a given system. In a work inspired by Gilfan ov et a l.| (^qI), 

into principal components. The resulting solutions are a posteri- ICervifio & Luridianal (|2006), demonstrate that EPS provides a 

ori constrained by additional assumptions to obtain the descrip- characterization of the probability distributions that describe all 

tion of the entire stellar population. In the second method (evo- the possible integrated luminosities of an ensemble of stars with 

lutionary stellar population synthesis, hereafter EPS ; see Tinsleil given evolutionary condition^ In general, this characterization 

[19681 for an example) the stellar populations are modeled fol- expressed in terms of the mean, p\{e,), of the distribution of 

lowing a prion rules related to stellar evolution and the relative t^e integrated luminosity of a system with a reference number 

contribution of each stellar population. Model results are com- ^jai-s, Nl,„^, that allows us to obtain the mean value of a dis- 

pared with observations, and properties of the stellar populations tribution for a different number of stars, N% by simple transfor- 

(and their evolutionary status) of the targeted source are inferred, mationfl p\ {U^-a, N*) = M). However, it is also possi- 

In both types of methods, degenerate solutions are present. The , .. ™ , , . ^ u j- . u 

, , , . , , „ , . , , ble to provide a more accurate characterization of such distribu- 

EPS method, however, IS based on stellar evolution theory and .• . , . r »u j- » u u »u 

. , , ', , , , , . , • , , tions using additional parameters of the distribution, such as the 
provides both a tool for both studying the integrated light of 

galaxies and comparing its properties with stellar model predic- i Throughout the paper, evolutionary conditions mean the age, metal- 

felons. licity, and star formation history of the stellar ensemble. 

^ In the following, we assume N„orm = 1 and we do not include the 

Send offprint requests to: M. Cervifio, V. Luridiana, L. Jamet; e-mail: subindex A in the luminosities to simplify the notation. We also omit the 

incs@iaa.es, vale@iaa.es, ljainet@astroscu.unam.mx subindex /i in the following equations. 
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second raw moment, ^^{(a), the variance, K2({a), the skewness, 

J 'l, and the k urtosis, 72, or even provide the distribution itself 
Cervino & L uridiana 2006). This natural description of EPS in 
probabilistic terms, in which it is possible to establish the dis- 
tribution of the integrated luminosity of a stellar ensemble, but 
not the precise position of this ensemble in the distribution, is 
not commonly considered in EPS literature, nor its application. 
A deterministic paradigm is usually applied: fi'^(nj)/N* is not a 
mean value but the actual fraction of stars of a given class present 
in the stellar system, regardless of the number of stars in the sys- 
tem. In the deterministic paradigm, given the total number of 
stars in the system and its evolutionary conditions, there is one 
and only one possible value of the integrated luminosity. 

However, the probabilistic description of stellar popula- 
tions, although not recognized and studied in detail before, 
has been present in the literature since the late eighties. 
iTonrv & SchneideJ (1198 8) presented the first description of the 
distribution of the integrated luminosities of stellar popula- 
tions making use of the mean and the variance in the asso- 
ciated theoretical distributions, and defined theoretical Surface 
Brightness Fluctuations (SBF, Ltot(A^*)) as the variance of the in- 
tegrated luminosity of a stellar population divided by its mean: 
^tot(A^*) = K2[Ltot{N*)]/fi'i[Ltot(N')]. Interestingly the depen- 
dence of K2[Lto((N*)] and |u'-^[L(ot(N*)] on the total number of 
stars cancels out numerically and Ltot(A^*) = ^ is independent of 
the number of the stars in the system. Assuming a Poisson dis- 
tribution for the number of stars in each given stellar population, 
SBF can also be expressed as the ratio of the second to the first 
raw moments of the distribution of integrated luminosities, and 
in this case, it can be redefined to be the luminosity-weighted 
lu minosity of the p opulation. 

iBuzzonil (Il989h presented an independent probabilistic de- 
scription of EPS, also based on the variance of the corresponding 
distributions, and defined the effective number of stars, N{€), to 
be the ratio ^\{(flK2{C). The value of 7V[Ltot(A^*)] for a system 
with A^* stars scales directly with N{() and it is an useful param- 
eter for estimating sampling effects produced by the discreteness 
of real stellar populations. In his formulation, it was assumed 
that the relative number of stars in a given population fluctuates 
around the th eoretical value way describe d by a Poisson distribu- 
tion (but see ICerviilo & Luridianall2006l and below). The con - 
nection between N and SBF was defined by iBuzzonil d 19931) . 
who related both quantities to the statistical entropy of stellar 
systems (see also Buzzoni 2008, for a recent review). 

As shown bv lTonrv & SchneideJ (1 19881) . SBF is not only a 
theoretical concept, but a real observable that provides a rela- 
tively direct technique for determining extragalactic distances. 
The amplitude of the relative fluctuations in the observed flux 
with respect to a suitable mean is related directly to the the- 
oretical SBF luminosity, and therefore to distance. In the ob- 
servational domain, its principal asset is being able to disen- 
tangle the fluctuations that provide physical information about 
the system from the noise. Since they are related to evolution- 
ary conditions, SBF provide informatio n about stellar evolu- 
tion and stellar population studies ( e.g. ICantiello et al. I l2003t 
iGon zalez-Lo pezlrra & Buzzonill2007l) . 

Most work performed on theoretical SBF was related to its 
calibration as a distance indicator (e.g. iTonrv et al. |[l990 ). al- 
though current calibrations are performed empiricall y. Further 



study exists on theoretical SBF themselves and their implication 
in the use of EPS: the very observation of SBF implies that EPS 
models cannot be interpreted in a deterministic way, but as a de- 
scription of a probability distribution that defines naturally an 
intrinsic metric of fit when physical properties are inferred frorn 
comparisons with model results (see lCervino & Luridianal2007l 
2008, for more details). 

On the other hand, the distributed nature described by 
EPS models was established clearly in theoretical studies of 
systems with a low number of stars, where the most prob- 
able number of stars in a given evolutionary stage «,/A^* 
differs significantly from the expected value i.e. the mean 
jUj(«,)/A^*. The characterization of integrated light in prob- 
abilistic terms began more r ecently. Partial solutions to the 
problems we re present ed by 'Cervino, Luridiana, & Castanded 
(2000); Bruzu ail (l2002h: iGirardi (2002); Lancon & Mouhcing 
(2000); Cervino etal. (2002); Cervino et al. (2001^ 
Cervino & Valls-Gabaud (2003); Cervino & Luridiana (.20041) . 
and Fagiolini. Raimondo. & Deg lTnnocentil ( l2007l) among oth- 
ers, and a comprehensive solution to the problem was presented 
by Cervino & Luridiana ( 2006), who introduced a probabilistic 
paradigm of EPS that applies to both under-sampled and well 
sampled systems. 

In this paper, we establish a framework for theoretical SBF 
based on the probability distribution theory of stellar popula- 
tions. The formalism allow us to describe the behavior of SBF 
for systems with a low number of sta rs per pixel, which was 
demonstrated by Aihar & Tonrvl (119941) to be a limitation of the 
method. Our main objective is to estabUsh a robust definition 
of SBF. In particular, we show that the Poisson assumption is 
not essential to the SBF definition. We examine the implications 
of this robust definition and explore new ways to ob tain obser- 
vational SBF. Paraphrasing Kernighan & Pikd ( Il999t) . "once the 
definitions are laid out, the theory tends to fall into place, and 
understanding it is comparatively easy". 

The present paper is organized as follows. In Sect. |2l we 
summarize the theoretical probability distributions involved in 
SBF studies and provide an introduction to the probabilistic de- 
scription of stellar populations. In Sect. [3] we present a detailed 
and quantitative description of SBF, establish the origin of SBF 
in probabilistic terms, and provide robust definitions of SBF ac- 
cording to the distributions involved. A discussion about the im- 
plications of the present result for SBF studies and additional 
SBF analysis methods are presented in Sect. |4] Finally, our con- 
clusions are given in Sect. |5] In a companion paper (Cervino et 
al. in preparation, hereafter Paper II), we will theoretically evalu- 
ate the observational error budget of the different ways to obtain 
SBF 



2. The distributions involved in SBF studies 



examp l es 01 aij ± ^ m s t ellar populations studies are IWorthe^ 
fm4): 'Condefli (1997); Br ocato et al. I (IT998h ; iBlakeslee et al. 
(2001 ), and Liu et al. ( 2000). In flie SBF definition, it is custom- 
ary to assume that the number of stars in a given stellar popula- 
tion follows a Poisson distribution. Unfortunately, no complete 



Following lTonrv & Schneid^d 19881) . observational SBF are ob- 
tained from a statistical analysis of fluctuations in the luminos- 
ity values of the pixels in the CCD image of a galaxy. First, a 
local mean is estimated at each point, and the fluctuations with 
respect to this mean are computed at each pixel. Second, the 
ratio between these fluctuations and the square root of the es- 
timated local mean is computed at each point. Third, the vari- 
ance in the resulting distr ibution is the SBF of the image. The 
ITonrv & Schneidedl 19881) method for obtaining SBF requires 
three distributions, although we will see in Sect. 14.21 that SBF 
can be obtained using only two distributions. We note that ^mc- 
tuation in Surface Brightness Fluctuations refers to the variance 
in renormalized fluctuations, i.e. the latter divided by the square 
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root of the local mean. This renormalization is needed to ensure 
that all pixels are sampling distributions with a common mean 
and variance. 

Although the term Surface Brightness Fluctuations is under- 
stood to be a measure of the moments of the luminosity func- 
tion from the analysis of the pixels in CCD images, we use a 
more general interpretation of SBF that relates any possible ob- 
servation of renormalized fluctuations to the stellar population 
properties of the system. 

The basic input is a set of observations of the same physical 
system. It can be the flux in a set of pixels from a galaxy (as 
proposed in the original method), a set of IFU obervations of a 
galaxy, a set of individual stars in a cluster or the integrated lumi- 
nosities of different portions of a system. The only requirement 
is that all elements sample either the same luminosity distribu- 
tion or a family of luminosities distributions whose parameters 
can be related to each other linearly (see below). In the follow- 
ing, we define resolution element to be any element of this sam- 
ple independently of whether we refer to, for example, a pixel, 
an IFU, an individual star 

The first distribution involved describes the possible inte- 
grated luminosities of a system for a given set of physical con- 
ditions, i.e. the luminosity distribution function of the popula- 
tion, or pLDF ((f(Lspj\f)). This distribution is directly related to 
the stellar luminosity distrib ution func tiorQI (sLDF: (p(()), by A^* 
successive convolutions (see lCervifio & Luridiana.2006> for de- 
tails). Implicitly, it is assumed that each resolution element is the 
result of a single random realization of ^(Lsp;A'*)- 

The second distribution involved is a reference distribution 
of local means (or, equivalently, expected values), needed to de- 
fine fluctuations. This is the key point of observational SBF stud- 
ies: when a system is observed, the number and class of stars in 
each resolution element are not a priori known, nor is the mean 
integrated luminosity corresponding to each resolution element. 
Hence, we need to estimate from the observation itself the ex- 
pected flux that would correspond to each resolution element as 
opposed to the actual flux of the element. This estimation can be 
performed only from the statistical analysis of a set of resolution 
elements. The expected flux estimate itself depends on arbitrary 
observational/analysis choices, such as the number of resolution 
elements used in the estimation, and the choice of the region 
used to obtain the expected value, and might vary from one ob- 
servation to another, or a data analysis to another (c.f. Sec. |4.2| l. 

Once the expected flux in each resolution element has been 
subtracted, the result is a pure distribution of fluctuations, that, 
by construction, has zero mean. However, the fluctuation of each 
resolution element depends on the number of stars in the element 
itself, which varies from element to elemenfl So it is needed to 
renormalize the fluctuation by the square root of the expected 
flux. By this operation, the resolution elements sample a set of 
distributions that, by construction, share the same mean and vari- 
ance. Observational SBF are just the estimation of the variance 
in this distribution of renormalized fluctuations. 

In this work, we differentiate between stellar population 
SBF, theoretical SBF, and observational SBF. Stellar population 
SBF are described in terms of the basic physical process respon- 
sible for intrinsic luminosity fluctuations across resolution ele- 
ments at fixed number of stars per resolution element. That is. 



This function provides the probability that a single star has a given 
luminosity, and it is only defined by the evolutionary conditions of a 
stellar sample. 

The higher the number of stars, the larger the possible fluctuation, 
as we see below (cf. Eg.ll3t. 
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the expected (mean) integrated flux in the resolution element is 
not obtained by observations but by a theoretical stellar popula- 
tion. 

Theoretical SBF are the SBF that can be theoretically pre- 
dicted for a given system. They can be shown to contain a con- 
tribution from stellar population SBF plus a contribution from 
the distribution of the number of stars per resolution element. 
That is, the expected integrated flux in the resolution element 
is obtained from a theoretical model of the distribution of ex- 
pected fluxes throughout the system. By their own definition, 
theoretical SBF measure a different physical quantity than stel- 
lar population SBF since they include an additional component. 
Finally, observational SBF are the observational counterpart of 
theoretical SBF. 

As anticipated in the introduction, our model of the physical 
origins of SBF is based on the probabilistic description of stellar 
populations. Some basic concepts of this field that will be used 
in the pape r are introduced in the follo wing subsection, which 
is based on lCervino & Luridianal (l2006l) and may be skipped by 
those readers who are already familiar with that paper 

2. 1 . The probabilistic description of stellar populations 

We introduce a few concepts used in this work. The probabilistic 
description of SBF is based on the concept of stellar luminos- 
ity distribution function sLDF, which is denoted by (pit), where 
i is the luminosity of individual stars. The sLDF provides the 
probability distribution of the luminosity of a randomly-selected 
single star. When referring to the sLDF, we will use /i,(/') to in- 
dicate central moments (or moments about the mean) and ;uj(/') 
to indicate raw moments (moments about zero). 

The explicit expressions for raw and central moments of the 
sLDF are: 



M = {"(fiOM^Vl'lpi, (1) 
Jo Y 

n„({) = ({-i^[rmde, (2) 

Jo 

where the mean jj\ {{) is the first raw moment, and the variance 
/i2(0 is the second central moment of the sLDF In these equa- 
tions, we have also expressed the integrals as a sum (as most 
synthesis codes do), where /, is the average luminosity of a given 
stellar type, the probability density that a randomly-selected 
star belongs to stellar type /, and A^st is the number of stellar 
types, such that 2^^^ p; - 1 . We note that the stellar types should 
be defined in such a way that the luminosity can be assumed con- 
stant or to have a small variance for each given stellar type, i.e. 
/,■ is a truly representative value of the particular stellar type; the 
traditional classification of stars into evolutionary stages, such 
as the main sequence, and giant branch, is invalid in this con- 
text because the stellar luminosity spans a wide range at each of 
these stages. 

We also use the cumulants k„ of different distributions. The 
fundamental aspects of the use of cumulant s in stellar population 
synthe sis were described extensively in Cerviiio & Luridianal 
(120061) and will not be repeat in detail here. Interested readers 
may find an introduction to this topic in any advanced textbook 
of probability and statistics (e.g. Kendall & Stuart 1977). We re- 
call only the relations between the first four cumulants and mo- 
ments: 
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and the definition of the skewness, y\, and the kurtosis, j2: 



71 



3/2 
A'2 

_ 



K3 

- — 
„2" 



(3) 
(4) 
(5) 
(6) 



(7) 
(8) 



When these symbols are expressed without a variate in 
parentheses or with the variate i, they refer to the sLDF. 

The same parameters can also be computed for the distribu- 
tion of the integrated luminosity of a population, pLDF, which 
contains A^* stars, (/j(Lsp;A'*): in this case, we write the depen- 
dence explicitly, e.g. //j(Lsp;A'')- Integrated luminosities are in- 
dicated in upper case, e.g. Lsp;^* is the luminosity emitted by 
A^* stars. Finally, we write the relation between the n-th order 
cumulant of the luminosity distribution of a cluster of A^* stars, 
KniLsp-N'), and the cumulants of the sLDF, k„: 



Kn(LsP-N') ^ N* X K„ 



(9) 



We make extensive use of this equation, which can be found 
in advanced statistical books o r, in the c ase of synthesis mod- 
els of integrated populations, in lCervino & Luridianal(l2006l) . We 
note that the cumulants of pLDFs (but neither the raw moments 
nor central nth-moments with n > 3) scale linearly with the num- 
ber of stars. In this work, we use K2 instead of fi2 to represent the 
variance of the distribution. 

The skewness, Fi, and kurtosis, F2, of the pLDF with A^* 
stars can be derived from Eqs.[3]-|9]and are given by: 



ri(isP;A'') 



Va^ 



7u 



T2(LsP;N-) = —72- 



(10) 
(11) 



Since for Gaussian distributions Fi = r2 = 0, Eqs.fTOlandfTTlin- 
dicate that the higher the number of stars A^*, the more Gaussian- 
like the distribution becomes. 

Finally, we use A^* to denote the number of stars per reso- 
lution element when this is a constant (stellar population case); 
N-^^ to represent the mean number of stars per resolution element 
in the area S used to estimate the mean luminosity in a resolu- 
tion element; A^pois to represent A^j;^, when the number of stars 
per resolution element is found to be represented by a Poisson 
distribution in the area considered; and A^ to refer to the number 
of stars per resolution element when this is a variable quantity. 

3. Theoretical definitions of SBF 

In Sect. 12] we emphasiz e the necessity of describ ing the theo- 
retical SBF definition of lTonrv & Schneideii(ll988l) arising from 
two factors. The first is the intrinsic dispersion in the total lumi- 
nosity of a population of a given number of stars A^; the second 
is the distribution of number of stars in our resolution elements. 



which is needed to define a reference value to define the fluctua- 
tion. 

The first of these factors is addressed in the theoretical 
framework of probabilistic synthesis models: a particular value 
of the integrated light is just a realization of all possible random 
combinations of a given number of individual stars with the same 
evolutionary conditions. Different realizations consist of diff'er- 
ent stellar mixefl and display different integrated luminosities. 
Even in the case of a strictly constant total number of stars in 
all resolution elements considered, fluctuations in the stellar mix 
always introduce a dispersion in the integrated luminosity val- 
ues. In the cases of an observation of an extended system with 
homogeneous evolutionary conditions, each resolution element 
provides the integrated luminosity of diff'erent realizations of the 
same (theoretical) stellar population. In our framework, we de- 
fine stellar population SBF to be those arising from this effect. 
The expression for population SBF is derived in the following. 

3.1. Population SBF 

We assume that all resolution elements have exactly the same 
number of stars, A^*. In such a case, the luminosity of each el- 
ement is a variable randomly-selected from the pLDF The first 
step in computing SBF consists of the characterization of the 
pLDF, ^(Lsp;Af)i the parent distribution of the L^p-^n', which is 
the luminosity of a resolution element with A^* stars. Equation|9] 
indicates that the mean luminosity and variance of the pLDF is 
related with the mean, yUj, and variance, K2, of the sLDF by: 

ti\(LsP:N') = A^Vi, (12) 

K2(L5P:N') = N*K2. (13) 

As we have emphasized before, the higher the value of A^*, 
the larger K2(Lsp-n-), or equivalently, the larger the amplitude 
of the possible fluctuations. Additionally, the theoretical scatter, 
cr(Lsp.^.) = s/ K2iLsp-,N') scales with yfN* and the mean inte- 
grated luminosity scales with A^*. Hence, it is natural to renor- 
malize the resulting fluctuation by -^///^(Lsp;^*) to eliminate nu- 
merically this dependence on A^*. The distribution of the renor- 
malized luminosity fluctuations, which is indicated by ^(Lg^^, ), 
is the distribution of the variable: 



LsP;N' - (Lsp-N' ) isP;A" - N* fl\ 



■^SP;N' 



^J^J.[iLsp■N•) 



(14) 



i.e. it can be obtained from (^(Lsp;A'*) by means of a linear trans 
formation. Using Eq.|4] the population SBF associated with reso 
lution elements containing A^* stars, i.e. the variance of ip{L^'^'.^, 
is: 



flue -J 



r ,..fluc ^ KliLsPM.) N*K2 



l<2_ 



(15) 



where Lsp;iv. is numerically independent of the number of stars 
per resolution element (assuming that all resolution elements in 
this region have exactly the same number of stars). Hence, we 



^ This is a simplification of the original theory, in which IMF sam- 
pling is but one of several ca uses of dispersion; it is , howe ver, sufficient 
to develop our argument. See lCervifio & Luridiand j2006l) for more de- 
tails. 
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can omit the subindex A^* and assume that the general expression 
for population SBF is: 



Lsp = -7. 



(16) 



However, although Lsp is independent of A^*, this does not 
mean that the shape of the distribution of renormalized fluctua- 
tions does not depend on A^*. The shape of this distribution can 
be inferred by computing its skewness and kurtosis. The 3rd and 
4th cumulants of (^(Lg™^.) are: 



N* K2 



K4 



and the skewness and kurtosis are therefore : 



/2 ' 



(17) 
(18) 



number of stars per resolution element, A^, in this region, i^s^ (A^). 
This distribution indicates the probability of having a resolution 
element with A^ stars, or, in frequency terms, the relative number 
of resolution elements with A^ stars when the number of resolu- 
tion elements is infinity. We assume that J ip^.^^ (N) dN = 1 . In the 
remainder of Sect. l3.2l we determine how the mean and variance 
of the fluctuation image depend locally on t//z^{N). These results 
are used in Sect. l3.3l to derive a general expression for SBF. 

The distribution of the possible luminosity values for the set 
of resolution elements considered, (f)(Lz^), is the result of the sum 
of all possible (^(Lsp;^) distributions, weighted by the distribu- 
tion of the number of stars A^ in the set of resolution elements 
defined by this area: 



(21) 



After simple but cumbersome algebraic operations, the fol- 
lowing is obtained: 



ri(^'sp'iv) 



r2(isp'Jiv) 



/C3 



4 



K4 



/3 \M 



K2 



-3/2 



VaF 



K2 



-2 



7 7'2- 



(19) 



(20) 



We note that Eqs. [19] and |20] provide the same values as 
Eqs. fTOlandfTTl so the shape in terms of the skewness and kur- 
tosis of y?(Lgp':^. ) is identical to that of the pLDF of a popula- 
tion of A^* stars, i^L,o,(i'SP;A'*)- Equations. [T9l and |20l also indi- 
cate that the shape of the distribution of renormalized fluctu- 
ations depends on the number of stars per resolution element, 
although all distributions share the same mean and variance. As 
an example, ICervifio & Luridianal ( l2006l) demonstrated that A^* 
must have values between 10^ (in visible bands) and 10'' (in in- 
frared bands) for Gaussian distributions of both the pLDF and 
renormalized fluctuations, to be achievable. That is, a given sys- 
tem would have a Gaussian distribution of renormalized fluctua- 
tions at one wavelength and a non-Gaussian distribution at other 
wavelengths. For a more detailed study of the number of stars 
needed to reach a Gausssiari re gime, we refer to the more de- 
tailed treatment presented in lCerviilo & Luridianal (l2006l) . 



3.2. Theoretical SBF: effect of a variable number of stars per 
resolution element in the estimation of local means 

The population Lsp (Eq. [TsT i is the theoretical SBF in the case 
in which all resolution elements have the same number of stars. 
In the general case, however, the number of stars per resolution 
element will differ for different elements and each element will 
sample a different pLDF. Hence, the local mean estimated from 
resolution elements in a region depends on the distribution 
of the number of stars in the resolution elements included in the 
region. 

We now derive the theoretical expression of SBF for the gen- 
eral case of estimation of a mean for resolution elements with a 
variable number of stars. In doing so, we must of course mention 
the set of resolution elements for which the local mean is com- 
puted, since this defines the distribution of the number of stars. 
We therefore indicate this region with E^, and introduce an addi- 
tional distribution that represents the distribution of the variable 



A<',(L2„) = fi\{N)fi\({), (22) 

^2(iE„) = yu'i (A^) ^2(^) + '<2(N) (23) 

XiiL^,,) = /"'i (A^) ^3(0 + 3^2(A^) yu'i (£) K2{() + K3(N) ^[\{), (24) 

K4(L^) = ^[(N) K4(e) + K2(N) [34(£) + 4i^[({) K^{i)] + 



6K,(N)i^[\()K2m + K4(N)fi[U£), 



(25) 



where a:,(A^) are the cumulants of the distribution i/'e^/A^), and 
Ki{{) are the cumulants of the sLDF 
The transformation 



(26) 



described previously now provides the distribution (p{L^^) with 
ti'AL^') = and /f2(4"') (i.e., the 'local' SBF) given by: 



K2{L.) 



and a skewness and kurtosis given by 



ri(L«-) 



r2(Lf^) 



K4iL^„) 
K2{L.)^' 



(28) 



(29) 



In general, the contribution from the variation in the num- 
ber of stars per resolution element in the region S^, in which the 
reference mean is deriveed must be modeled and evaluated for 
each system if one wants to compare the SBF with stellar popu- 
lation SBF. However, in the case of relatively close galaxies with 
smooth luminosity profiles, it can be assumed that the luminos- 
ity profile in any region is almost flat and that ^•e,^{N) follows 
a Poisson-like distribution. We consider this specific case in the 
following section. 
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3.2.1 . Case of Poisson distribution in tlie number of stars 

There are two cases in which the above equations provide a sim- 
ple relation to the moments of the sLDF: (a) ifJi-^iN) is a Dirac 
delta function centered on A^*; (b) i/zz^X^) is ^ Poisson distribu 
tion with mean value A^, 
to the A? = N* 

case, it is known that fJ.\{N) = K2{N) = Ki{N) - K4{N) = N, 
Introducing these relations in Eqs.|22]to|25] we obtain: 



pois- The first of these cases corresponds 
constant case discussed above. In the Poisson 



pois ■ 



/^'i(^2.;iVpoJ = ^pois/j'i, (30) 

K2{L^,;Nf,-J = A^poisfe +/^f ) = A^poisjU2' (31) 

KiiL^^rJf^J = A^pois/^s, (32) 

f4(i'I„;A'p„J = A^pois/i4- (33) 

These equations characterize the distribution of luminosi- 
ties in the case in which the distribution of the number of stars 
per resolution element follows a Poisson distribution with mean 
A^pois- We note that K2{Lj_^-^t^^.^) contains a component equal to the 
variance of a pLDF with ^pois stars, plus a component equal to 
the variance in luminosity of a Poisson distribution of identical 
stars with individual luminosity given by ji'^ . We can characterize 
the distribution of renormalized fluctuations <p(L^\, ) by: 



^1 

r /rfluc s _ 1 ^3 



'/ii'^pois' AT ,,'2 

-'»pois 



(34) 
(35) 
(36) 



Equation [34] is a rig orous proof of the assertion of 
iTonrv & Schneided ( Il988h that the SBF can be calculated to be 
the ratio of the first to the second raw moment of the stellar 
luminosity function for a Poissonian approximation to the dis- 
tribution of the number of stars per resolution element used to 
compute the mean flux of each resolution element. However, we 
note that a Poissonian distribution of the number of stars per res- 
olution element must not be taken for granted in all situations. 

3.3. Theoretical SBF: global effect of a variable number of 
stars in the resolution elements 

In the previous case, we obtained the distribution of possible 
renormalized fluctuations in a 'local' region that includes a given 
resolution element and where the mean had been obtained from 
the analysis of an additional set of resolution elements defining 
a region S^,. In this section, we derive the expression for theoret- 
ical SBF by assuming that the local luminosity mean values are 
obtained as above and integrating the local SBF over the entire 
system. Generally speaking, the distribution of stars per reso- 
lution element in each E^, throughout the galaxy will vary and 
there will be a different set of Eqs. |22] to |25] for each resolu- 
tion element. That is, we have a i//-i:^{N) distribution defined for 
each resolution element. Since the problem is difficult to solve in 
its generality, to simplify we assume that the four cumulants of 
^"LfiiN) can be expressed as linear functions of the mean number 
of stars per resolution element, A^j^, inside each 2^ (which may 
vary from element to element). Implicit in this hypothesis is that 
the cumulants that depend on the stellar population (yUi(^), K2i(), 



etc.) are the same for all resolution elements, i.e. the stellar pop- 
ulation is homogeneous across the region in which the final SBF 
are computed, or the sLDF is universal, or any variation in the 
stellar population properties occurs smoothly. We note that this 
hypothesis is implicit in any SBF analysis and is not a particular 
assumption of this work. 

With this hypothesis, we define 2sbf to be the region in 
which final SBF are evaluated; it can be the entire or a region 
of the galaxy. The distribution of renormalized fluctuations in 
which we are interested, <l)2j,gp(L'''"^), is the result of the sum of 
all possible 0(L^'"^) distributions weighted by the distribution of 

N^^, given by *I'sbf(A^i;^), which is related to the luminosity pro- 
file of the galaxy region over which the SBF is computed. 

The distribution of renormalized fluctuations can be obtained 
as in the previous case: 



*I'sBF(A^L,)0(L^:^)fl'A^E, 



(37) 



We note that this equation is formally equal to Eq. [JT] but in 
place of the local distribution of the number of stars from Eq.l2n 
we have the global distribution of the mean number of stars. 
Furthermore, in Eq.|2T]we had the pLDF, whereas in Eq.[37]we 
have the distribution of renormalized fluctuations 0(L^"'^). This 
is necessary because of the limitations intrinsic to the observa- 
tional method: since the local mean and the distribution of stars 
per resolution element are unknown, the local mean must be ob- 
tained from a set of resolution elements, then used to obtain the 
fluctuation image (which depends on the definition of the local 
area S^), with which, finally, the SBF are obtained. 



Case of constant number of stars in all resolution elements. 
Let us obtain the cumulants of this distribution in the theo- 
retical case of a fixed number of stars per resolution element 
4>{L^'^) - l/;L,„,(^'sp;^,.)■ We must perform the same algebraic op- 
erations as in the previous case, but by using instead the fluctu- 
ation distribution of the pLDF In the case of A:2[<l'j;sBp(Lgp^,)], 
we find: 



i = ^2[1>s(isP;iV.)] = =isP, 
^1 

and the skewness and kurtosis of <l>(Lj™^,) are: 



(38) 



ri[<i)E(4p;iv.)] = n 



SBF 



(39) 
(40) 



SBF 



where the bracket notation < x >sbf denotes the expected value 
of X over the region in which the SBF is computed. 

Case of Poisson-distributed number of stars. If we assume 
that the number of stars per resolution element in the regions 
we use to obtain the mean values follows a Poisson distribution, 
(piL^^) — (f){L^\, ), the following relations are obtained: 



flue 

'pois;A'p, 



_ ^"2 _ - 



(41) 
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A'3 / 1 



'pois;A'pois-'-l /3/2 



1^2 ^ y/^pois I 
/ 1 



\"P™s/sBF 



(42) 
(43) 



In both cases, the shape of <t>j:(L^"'^) depends on <1/Va^) 
and (l/N). We note that these values are in general differ- 



ent from l/-\/(N) and 1/(N). In particular, it can be shown 
jKendall & Stuartill977l exercise 9.13) that the product <A^) x 
(1 /N) > 1 if > 0, as in our case. The equality only holds if the 
distribution of is a Dirac delta function. Hence, in general, the 
use of 1/(A^) provides only lower limits to the estimation of the 
shape of ^(L*'""). 



Implications in the general case. An analogous reasoning, but 
leading to more complicated expressions, applies to the general 
case. There are two relevant aspects here: (1) The global SBF 
coincides with the local SBF. Since the latter contains the contri- 
bution of the stellar population SBF plus the contribution of the 
variation in the total number of stars in the region used to obtain 
the reference mean value, the final SBF also includes these con- 
tributions. (2) In general, the distribution of renormalized fluc- 
tuations, ^^(L'^"'^), does not follow a Gaussian distribution and 
its shape will depend on the value of < 1 / > and < l/N > 
and the raw moments or cumulants of the sLDF Only for a suf- 
ficiently high number of stars per resolution element will this 
distribution be Gaussian. A qua ntitative evaluation of how high 
this number can be, is found in Cervifio & Luridianal ( 1200 6'). As 
reference values, Fi and Fa must be lower than 0.05 to achieve 
a Gaussian-like symmetric distribution (at a 3cr level), and Fi 
and F2 must be lower than 0.5 to avoid a bimodal distribution 
of renormalized fluctuations. As an example, a Gaussian distri- 
bution of renormalized fluctuations is obtained for systems with 
more than 10^ stars per resolution element (e.g. pixel) for old 
stellar clusters in infrared bands. 

We now explore the shape of the distribution of renormalized 
fluctuations for the z band used in HST-ACS surveys (Mei et al. 
l2007h . Assuming that the typical values of skewness and kur- 
tosis of the sLDF for the z band of an old stellar popula- 
tion are simillar to those o f the I band, we use the results of 
ICervifio & Luridianal (I2006h where ri(/) ^ 10^ and joil) « 10\ 

For the Virgo cluster observed by the ACS/WFC, the pro- 
jected distance in a pixel is » 4 pc assuming a distance of D » 
16.5 Mpc and a pixel scale x 0.05"/pix. Assuming an aver- 
age p rojected density of stars of 10^* stars pc~^ jTuntsov & LewisI 
I2OO6 ). the number of stars per pixel is (A^) x 1 .6 x 1 0^. This num- 
ber yields values of Fi[(t>j:(/)] ^ 0.8 and Vildi^d)] ~ 0.06. That 
is, the distribution of renormalized fluctuations is asymmetric 
and possibly multimodal. It is also possible to use a larger resolu- 
tion element for SBF computation corresponding more closely to 
the ACS/WFC z band PSF, which can be approximated by areas 
of 3x3 pixels. In this case, both Fi[il)£] and F2[<1'2;] decrease by 
factors of 3 and 9 respectively, which yields Fi [<I>2;(/)] ~ 0.26, 
still a factor of 5 higher t han that required to ensure gaussianity. 

Considering M32 an d iTonrv & Schneider ( 1988), we assume 
a distance of 0.9 Mpc and a pixel scale of x 0.415"/pix; 
the projected distance in a pixel is then 1.8 pc, which yields 
(N) ^ 3.5 X lOl In flie V ban d, with yi(V) ~ 50 and y2{V) ~ 400 
( iCervifio & Luridianal I2OO6I) , it is found Fi[<Di:(y)] ~ 0.8 and 
r2[<l>5;(y)] ^ 0.1. This distribution is again non-Gaussian, in 
contrast to the assumption made by iTonrv & Schneidedl (Il988l) 
based on an estimated frequency of 20 giant stars per pixel. We 



note that giants are a broad stellar type, which cover a broad 
range in luminosity; taken at face value this number, does not 
therefore provide sufficient information about the shape of the 
distribution of renormalized fluctuations. 



4. Discussion 

Once the distribution involved in SBF has been obtained, it is 
possible to propose a robust definition of theoretical SBF, es- 
tablish guidelines for comparisons between theoretical and ob- 
servational SBF, propose additional methods to obtain SBF and 
extend its range of applications. In this section, we first discuss 
some issues related to SBF definitions found in the literature 
(in particular the unnecessary mention of Poisson distributions). 
Secondly, we propose additional methods to obtain SBF and dis- 
cuss the advantages and differences. Finally, we show some ad- 
ditional applications of the method. 

4. 1. Poisson or not Poisson? 

The most common explanation of the origin of SBF found in the 
literature is Poisson fluctuations in the number of stars (and in 
the number of stars per stellar type in a given population) across 
different resolution elements. However, we have shown that in 
this type of definition two different concepts are entangled: 

1 . The first concept involved is the population SBF, which ap- 
pear naturally as a manifestation of the distribution of the 
possible luminosity of any ensemble with a fixed number of 
stars and physical conditions. 

2. The second concept involved is the distribution of the to- 
tal number of stars in different resolution elements that are 
needed to obtain a reference value for the definition of fluc- 
tuations. We have shown that it can be understood as the pos- 
sible fluctuation when a population consists only of a partic- 
ular type of stars, all of the same luminosity (. That is, the 
stellar population itself does not produce any fluctuation at 
all, and the only component of the SBF arises as a conse- 
quence of variations in the total number of stars. In this case, 
the SBF would be { if the total number o f stars per resolution 
elements follow a Poisson distribution jTonrv & SchneideJ 
mH), although a more general, exact, and simple formula- 
tion can be obtained without the use of Poisson approxima- 
tion. 

In our work, we have assumed that (1) the number of stars 
per resolution element is fixed (and produce the stellar popula- 
tion SBF) although (2) it is impossible to determine the number 
of stars per resolutio n element and we must co nsider instead an 
average value. When lTonrv & Schneideij(ll988l) defined theoret- 
ical SBF on a Poisson base, they used these two different con- 
tributions implicitly: that of the stellar population SBF and that 
due to the local variation in the mean number of stars per resolu- 
tion element (where the Poisson distribution applies), but with- 
out a clear distinction between both components. We emphasize 
again that Tonrv & Schneider ( 1988) make use of a probabilistic 
description of EPS since they assume that the number of stars 
of a given stellar type does not have a fixed value, but fluctuates 
about a mean. Under a deterministic paradigm of EPS, where 
for given evolutionary conditions and total number of stars there 
is one and only one value of the integrated luminosity, the only 
possible source of fluctuation is that related to the total number 
of stars. This would explain the customary mention to Poisson 
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fluctuations in SBF definitions and its lower relative impact on 
the EPS literature. 

It is now necessary to evaluate the contribution of the 
Poisson component to the theoretical SBF; only in the case that 
it provides a measurable contribution, it is justified to mention it 
in the definition. 

To evaluate this contribution, we use N defined by Buzzoni 
(note that in its original definition N is defined on a Poisson 
basis, so it uses the second raw moment), and compare this with 
i-TSSS = i'j;^;A'„„i, defined in Eq.[34l 



1 

77 



A*, 



K2 + fJ-f 

~7r 



+ 1. 



Hence, the relationship between Lsp and Lts88 is 



^SP 



i'SP 



= 1 - AT. 



(44) 



(45) 



;pois 



N has values ar ound of 10"^ (iBuzzonil 1 19891 119931: 
ICervifio et al. 1 12002|) when the stellar luminosity distribution 
function is normalized to an initial total mass in the population 
of 1 Mq, and has a similar order of magnitude when normal- 
ized in terms of the number of stars. The numerical values of 
the pure population SBF, Lsp, and population SBF plus Poisson 
contribution due to the variation in the number of stars per res- 
olution element when the reference mean is obtained, Ltsss, are 
therefore similar. That is, Poisson fluctuations in the total num- 
ber of stars only contribute to around 0.1% of the total SBF. As 
a consequence, it is very ambiguous to make an explicit mention 
in the definition of SBF a contribution that only accounts for a 
0.1% of the variance, and not mention the intrinsic distributed 
nature of the integrated light of a stellar population that is re- 
sponsible for 99.9% of the SBF. Fortunately for previous SBF 
studies based on EPS, this change in the definition of SBF has 
a negligible numerical impact on current SBF computations of 
Single Stellar Populations models. 

This result also shows that (stellar populations) SBF should 
be an intrinsic result in EPS studies, and not a particular compu- 
tation designed for a restricted range of observational cases. 

4.2. Different ways to obtain SBF 

Once the nature of the different components implicit in theo- 
retical SBF has been established, it is possible to define different 
ways to compute SBF. The methods presented here differ mainly 
in the way the resolution element is defined and how the refer- 
ence mean is obtained. 

The first is that established by iTonrv & Schneiderf (Il988h . 
which has been the reference method used in this paper. Its main 
characteristics are: (a) the resolution elements are the individual 
pixels in a CCD image; (b) the reference mean is obtained from a 
smooth model image of the galaxy profile with a posterior t reat - 
ment of large-scale residuals (we refer to iMei et al1l2005L and 
references therein for further details and improvements); and (c) 
the SBF are obtained by means of the spatial power spectrum of 
the entire galaxy or a galaxy region (in general annular regions 
around the center). In this method, the main advantages are the 
use of a large number of resolution elements (the individual pix- 
els), and the use of the power spectrum analysis, which allows a 
clear discrimination between noise (with a flat power spectrum) 
and SBF signal. Additionally, the method includes other sources 
of fluctuations (such as the contribution of globular clusters) that 



we have not considered in this work. However, the application 
method is far from simple and there are several observational 
aspects that must be taken into consideration. Intereste d read 



ers may find more detail, for examples, in Tonrv & Schneidei 



('l988').'Aiha r & Tonrvl jl994 . lTonrv et alj tl 997,) . Jensen et al 
(2003), Mei e t all (120051) and references flierein. 



A second method is that proposed by Bu zzonil (l2008 b: (a) 
The image is divided into several regions following the symme- 
try of the system (i.e. pie-like sections) and the flux is obtained 
in each section. These sections (that may contain a large num- 
ber of pixels) are the resolution elements; (b) The average flux is 
obtained from all the resolution elements and the renormalized 
fluctuation with respect to this common mean is computed; (c) 
The variance of the resulting distribution of renormalized fluc- 
tuations is obtained. In this case, the resolution element clearly 
includes both the SP component and the profile component, al- 
though in fact it is just a natural extension of the original method 
for which, obviously, the Poisson approximation does not apply. 
Although this method may not provide accurate SBF values, it 
can be used as a first SBF estimation and for stellar populations 
studies (provided that the SP component is more important than 
the profile component). 

Other ways to obtain SBF can be propose d, such as 
t he application to resolved systems ( Aihar & To nryl [19941 ; 
iRaimondo et al. 1 120051) . the analysis of the ensemble of IFUs 
observations of the system as opposed to the analysis of each 
particular IFU, a combination of both methods, or, even, the use 
of data archives with different observations of the same system. 
The crucial points to remember in these situations are: (a) how 
the stellar population is sampled in each one of the considered 
resolution elements and (b) how the distribution of the mean 
number of stars in the resolution elements is used to define a 
mean reference flux for each element. 



4.3. Stellar population SBF as an intrinsic metric of fit in EPS 
studies 

In the previous section we provided different ways to obtain SBF 
from the data. However, we showed that SBF are an intrinsic re- 
sult of the modeling of stellar populations that provides informa- 
tion on t he dispersion in the integrated luminosity of any stellar 
system (iBuzzonilll993L 12008 ). Furthermore, stellar population 
SBF are a scale-free quantity independent of the number of stars 
in the system (although the shape of the distribution of renormal- 
ized fluctuations is not scale free). Hence, stellar population SBF 
provide an intrinsic metric offlt for any comparison of a single 
observation (e.g. the inte grated spectrum o f a galaxy) with syn- 
thesis models results (see lCerviflo & Lurid iana 20071 l2008l for 
more details). 

When a single observation is compared with the mean ob- 
tained by EPS models, it is possible to obtain its associated 
renormalized fluctuations precisely as achieved in Sect. 13.11 us- 
ing Eq. [14] Hence, we can obtain the position of the individ- 
ual observation in the distribution of renormalized fluctuations 
i^(Lgp^.) and establish an associated probability for the corre- 
sponding fit. In other words, stellar population SBF provide a 
tool to establish the accuracy of the fit, as opposed to current 
methods that only provide precise values without any accuracy 
information. 

In particular, the comparison of stellar population SBF at 
different wavelengths infers, that it is impossible to fit data at 
infrared wavelengths with the same precision as at optical wave- 
lengths, since the first have a higher SBF value, and hence a 
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larger intrinsic dispersion. This result is independent of the num- 
ber of stars in the observation. As a consequence, a probabilistic 
paradigm in the use of synthesis models is necessary not only for 
the analysis of systems with a low number of stars (affected by 
IMF sampling effects) but for an accurate analysis of any system, 
regardless of size. 



5. Conclusions 

In this work we have investigated the nature of SBF. We have 
introduced a distinction between stellar population SBF, theoret- 
ical SBF, and observational SBF. Stellar population SBF depend 
only on the intrinsic properties of stellar populations; theoretical 
SBF are determined by both the population and the contribution 
from the distribution of the total number of stars in the used res- 
olution element; and the observational SBF represent the obser- 
vational counterpart of theoretical SBF. Stellar population SBF 
and theoretical SBF should be considered in probabilistic terms. 
On the other hand, observed SBF are statistically distributed and 
should therefore be interpreted in statistical terms. 

In this paper, we have applied the probabilistic description 
of synthesis models to SBF. We have found that SBF are sim- 
ply a manifestation of the way in which Nature samples either 
the integrated luminosity distribution function (pLDF) in inte- 
grated systems across different resolution elements, or the stel- 
lar luminosity distribution function (sLDF) of resolved systems. 
We have also found that this definition does not require any as- 
sumption about the distribution of stars across different stellar 
types. Stellar population SBF are therefore simply defined to be 
the ratio between the variance and the mean of the luminosity 
distribution function in linear fluxes: 



_ K2 _ K2(Ltov,N') 

i-SP - — - —TTi \- 

Hi HiiLtot-N') 



(46) 



In addition to the definition of stellar population SBF, we 
have introduced a definition for theoretical SBF, which con- 
sists in general of a theoretical description of an observational 
method. We have shown that the effect of the distribution in 
the number of stars per resolution element, A^, is included in 
theoretical SBF in two ways. First, this distribution affects the 
computation of the mean flux needed to define the amplitude of 
renormalized fluctuations for each resolution element. Second, 
the distribution of mean fluxes also affects theoretical SBF. Both 
effects introduce a term that is added to the population contri- 
bution. Only for flat galaxy profiles can this term be a simple 
Poisson contribution; in the case of a different profile, the term 
differs and must be evaluated explicitly. As a result, theoretical 
SBF are defined to be: 



_K2 , K2{N) _ - , K2(N) 



(47) 



Observational SBF are an estimation of theoretical SBF and 
must be corrected for the contribution of the variation in the total 
number of stars per pixel, which depends on both the observa- 
tional strategy (by means of the way in which the local mean 
is computed) and the particular galaxy (by means of the galaxy 
profile), before being compared with stellar population SBF. 

Assuming, as a particular case, that the total number of stars 
per resolution element is Poisson-distributed, we have shown 
that the corresponding contribution is negligible compared with 
the contribution of the stellar population in the break-up of theo- 
retical SBF. Therefore, the common theoretical definition of SBF 



as "mean, luminosity-weighted luminosity of the stellar popula- 
tion" resulting from a Poisson assumption is a restrictive defini- 
tion that does not necessary apply to the observed SBF, nor any 
definitions based on an interpretation of SBF as Poisson fluctua- 
tions in the number of stars. 

We have also characterized the distribution of renormalized 
fluctuations for different cases in terms of skewness and kur- 
tosis, showing their relation with the skewness and kurtosis of 
the sLDF We have shown that the shape of the renormalized 
fluctuation distributions for a system with A^* stars per resolu- 
tion element is similar to the shape of the integrated luminosity 
distribution of clusters with A^* stars, the population luminosity 
distribution function (pLDF ). Taking advantage of the results by 
iCerviiio & Luridianal (l2006h . we have shown that the renormal- 
ized fluctuations of flux across galaxies are not Gaussian, as had 
been commonly assumed. 

Additionally, we have shown alternative ways of compute 
SBF and that stellar population SBF define an intrinsic metric of 
fit for the comparison of single observations with synthesis mod- 
els results. In fact, SBF are observationaling evidence that the 
results of evolutionary synthesis models cannot be interpreted in 
a deterministic way, but must be interpreted using a probabilistic 
paradigm. 

Finally, as a caution to theoretical computations of SBF, it 
is fundamental to ensure that the stellar types considered in the 
computations are defined in so that the luminosity can be as- 
sumed to be constant within each stellar class considered. SBF 
computed with isochron es that in clude a detailed treatment of 
the post-AGB evolution jCioni e t al. 2006a b) yield only lower 
limits to real SBF, since the post-AGB sequence provided by 
these isochrones are the average luminosity during the post- 
AGB evolution. The variance in this average luminosity must 
also be known before realistic SBF computations can be ob- 
tained based on these isochrones. 
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